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Abstract 

We present a detailed analysis of the spin models with near-neighbors interactions con- 

^^ structed in our previous paper [1[ by a suitable generalization of the exchange operator 

(^ . formalism. We provide a complete description of a certain flag of finite-dimensional spaces 

^^ I of spin functions preserved by the Hamiltonian of each model. By explicitly diagonalizing 

I the Hamiltonian in the latter spaces, we compute several infinite families of eigenfunc- 

Cin' tions of the above models in closed form in terms of generalized Laguerre and Jacobi 

■^ ■ polynomials. 
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1. INTRODUCTION 

The discovery of the quantum models named after Calogero [2] and Suther- 
land [3|, |j] is a key development in the theory of integrable systems which has ex- 
erted a far-reaching influence on many different areas of Mathematics and Physics. 
This is borne out by the relevance of these models in such disparate fields as 
groiip theor y Isl , |6|, the theor y o f special functions and orthogonal polynomi- 
als [71, Isl, y, |lO|, soliton theory [ll[, random matrix theory [l^, |l3|, |lj], quantum 
field theory [l5|, |l6|, |l7|, etc. The Calogero and Sutherland models describe a sys- 
tem of A^ quantum particles in a line or circle, respectively, with pairwise interactions 
inversely proportional to the square of the distance. Over the years, many different 
generalizations of these models have been considered in the literature. One such 
significant extension was proposed in the early eighties by Olshanetsky and Perelo- 
mov [l8i], who showed that both the Calogero and Sutherland models are limiting 
cases of a more general integrable model with a two-body interaction potential of 
elliptic type. The integrability of the latter model was explained by expressing the 
Hamiltonian as one of the radial components of the Laplace-Beltrami operator in a 
symmetric space associated with the A^r-i root system. It was also shown in Ref. [l8| 
that one can construct integrable generalizations of the Calogero-Sutherland (CS) 
models associated with any classical (extended) root system, like BCn- 

Another essential feature of the original Calogero and Sutherland models and 
their generalization to other root systems is their exact solvability, that is, the fact 
that the whole spectrum can be computed in closed form using algebraic techniques. 
In the last decade, some authors have introduced further extensions of CS models 
which are quasi-exactly solvable, in the sense that only part of the spectrum can 
be computed algebraically [l9|, |20|, |2l|, |2^. In all of these quasi-exactly solvable 



CS models, the Hamiltonian can be expressed as a polynomial in the generators of 
a realization of the Lie algebra si^+i in terms of first-order differential operators. 
Since these operators leave invariant a finite-dimensional space of functions, the 
Hamiltonian is guaranteed to possess a finite number of eigenfunctions belonging to 
this space. 

A great deal of attention has also been devoted to constructing models of CS type 
for particles with internal degrees of freedom (typically spin), partly motivated by 



their intimate connection with integrable spin chains of Haldane-Shastry type [23|, 
24| . Two main approaches have been followed in order to incorporate spin into CS 
models, based either on supersymmetrv 25l . l26l . l27l. l28l. l29l| or the exchange (also 
known as Dunkl) operator formalism JSJ, |30|, ISjJ, |32, |33|, IS^]- The spin models thus 
obtained include the exactly solvable spin counterparts of the scalar Calogero and 
Sutherland models of Aj^ and BCn type, as well as several quasi-exactly solvable 
deformations thereof, some of them with elliptic potentials |35l . |36|. A common 



property shared by all of these models is the long-range character of the interaction 
potential, in the sense that all particles interact with each other. 

The connection between spin CS models and spin chains of Haldane-Shastry 
type was first elucidated by Polychronakos through a mechanism known as the 
"freezing trick" [STj. The main idea is that in the large coupling constant limit 
the particles in a (dynamical) spin CS model freeze at the classical equilibrium of 
the scalar part of the potential, thus giving rise to a spin chain with long-range 
position-dependent interactions. In this limit the eigenfunctions of the spin CS 



model factorize into the product of an eigenfunction of the corresponding scalar CS 
model times an eigenfunction of the associated spin chain. If all the eigenf unctions 
of both the scalar and spin CS models are known, the partition function of the 
corresponding spin chain can be exactly computed from those of the scalar and spin 



CS models \m M 40 



A few years ago, Jain and Khare presented a novel class of scalar CS-like mod- 
els of A]sf type, characterized by the fact that each particle only interacts with its 



nearest and next-to-nearest neighbors [41|. In a subsequent paper [42|, Auberson, 
Jain and Khare discussed a generalization of these models to the BCn root system 
and to higher dimensions. The latter papers, however, left open some important 
issues, such as the exact or quasi-exact solvability of these models, the derivation of 
general explicit formulas for their eigenfunctions and the existence of similar mod- 
els for particles with spin. The last question was first addressed by Deguchi and 
Ghosh [43], who introduced and partially solved several spin 1/2 extensions of the 
scalar models of Jain and Khare using the supersymmetric approach. By a suitable 
generalization of the exchange operator formalism, in our previous paper [l] we con- 
structed the three spin models of Aj^j type with near neighbors interactions listed in 
Eq. ([2]) below. A significant property of these models is the fact that the spin chains 
obtained from them by the freezing trick feature short-range position-dependent in- 
teractions, and thus occupy an intermediate position between the Heisenberg chain 
(with short-range position-independent interactions) and the spin chains of Haldane- 
Shastry type (possessing long-range position-dependent interactions). In Ref. [l| we 
presented without proof closed- form expressions for several infinite families of eigen- 
functions of the scalar reductions of all three models, considerably generalizing the 
results of Ref. 42] . We were also able to derive similar expressions for a wide class 



of spin eigenfunctions of the models fl2b|) and fizcl) . The computation of the spin 
eigenfunctions for the remaining model (l2all . which is probably the most interest- 
ing one due to the rich structure of its finite-dimensional invariant spaces, was not 
undertaken in Ref. pb]- 

In this paper we present a detailed analysis of the models ([2]), with special empha- 
sis on the rational model (l2a]) . In particular, we have achieved a complete description 
of the flag of invariant finite-dimensional spaces for the latter model presented in 
Ref. [l|. More importantly, we have found that this flag can be further enlarged with 
an additional family of spin functions. We have computed all the eigenfunctions of 
the model (!2aj) belonging to the new flag, thereby obtaining seven infinite families 
of spin eigenfunctions in closed form. These eigenfunctions have been expressed in 
all cases in a compact way in terms of generalized Laguerre and Jacobi polynomials. 
The resulting expressions will be used in a forthcoming paper for computing a num- 
ber of eigenvalues and eigenfunctions of the spin chain obtained from the model (l2al) 
by taking the strong coupling limit. 

The paper is organized as follows. In Section [2] we define the Hamiltonians of 
the spin many-body models which are the subject of this work, and show that 
they can be expressed in terms of suitable differential operators with near-neighbors 
exchange terms. Section [3] is entirely devoted to the characterization of certain finite- 
dimensional spaces of polynomial spin functions invariant under these operators. 
The first part of this section deals with the construction of the latter spaces and the 
proof of their invariance, cf. Theorem[T]and CorollarylH In the rest of the section we 
complete the description of the invariant spaces for the model (!2al) . by identifying the 



spin states which satisfy a restriction stated in Theorem [H In Section H] we show that 
the eigenvalue problems for the Hamiltonians of the models ([2]) restricted to their 
invariant spaces reduce to finding the polynomial solutions of a corresponding system 
of differential equations. By completely solving the latter problem, we obtain several 
(infinite) families of eigenf unctions for the models ([2]), whose explicit expressions 
are presented in Theorems [2H11 Finally, in Section we summarize our results and 
outline some related open problems. 

2. THE MODELS 

In this section we shall introduce the three types of A^-body models with near- 
neighbors interactions whose study is the aim of this paper. We shall also recall 
from our previous paper [1|] the relation of each of these models with a corresponding 
differential-difference operator involving near-neighbors exchange operators. 

Let us begin with some preliminary definitions. We shall denote by |si . . . sn), 
where Si = —M, —M + 1, . . . , M and M G |N, the elements of the canonical basis 
of the space E of the particles' internal degrees of freedom (SU(2M + 1) spin). The 
action of the spin permutation operators Sij on this basis is given by 

Oij I . . . 5j . . . 5j . . . / I . . . 5j . . . 02 • • • / • 

The operators Sij can be expressed in terms of the fundamental SU(2M + 1) gener- 
ators S^,a = l,.. . ,4M(M + 1), as Sij = 1/(2M + 1) + ^„ ^fS";. The Hamiltonians 
of the models we shall be concerned with are given by 

(1) 
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a-Si^i+i), (2a) 



[a-Si^i+i], (2b) 

V2 = 2a^ ^ cot(xi - Xi-i) cot{xi - Xj+i) + 2a ^ csc^(xi - Xj+i) (a - Si^i+i), (2c) 

i i 

with r^ = "^ixf and a,b > 1/2. Here and in what follows, all sums and products 
run from 1 to A^ unless otherwise stated, with the identifications Xq = x^ and 
xat+i = Xi. a few remarks on the configuration spaces of these models are now in 
order. In all three models the potential diverges as (xj — Xj+i)"^ on the hyperplanes 
Xi = Xj+i, so that the particles i and i + 1 cannot overtake one another. Since we 
are interested in models with nearest and next-to-nearest neighbors interactions, we 
shall henceforth assume that xi < ■ • • < xn- For the second potential (l2bl) we shall 



take in addition xi > 0, due to the double pole at Xi = 0. For a similar reason, we 
shall assume that Xj+i — Xi < it for the potential (l2cl) . 



Remark 1. The Hamiltonians ([Tl) admit scalar reductions Hf^ ^ HA satisfy- 

ing the obvious identity 

H,{Ms)) = {HTns) . 

where ^/^ is a scalar function of the coordinates x = {xi, . . . ,xm) and \s) is a 
totally symmetric spin state. It follows that the spin Hamiltonians H^ possess 
factorized eigenf unctions of the form \1' = ^|,s), where ip is an eigenfunction of 
the corresponding scalar Hamiltonian H^'^ and \s) is again a symmetric spin state. 
The scalar reductions of the models ( !2a|) and ( !2cl) were introduced by Auberson, 



Jain and Khare in Ref. [42f|, whereas that of the model fl2b[) first appeared in our 
paper [l|. It should also be noted that for spin l/2_the potentials Vq and V2 differ 
from those studied by Deguchi and Ghosh in Ref. 1431] by a spin-dependent term. 

The models (|2]) share a common property that is ultimately responsible for their 
partial solvability, namely that each Hamiltonian H^ is related to a scalar differential- 
difference operator involving near-neighbors exchange operators. Indeed, let Kij 
denote the operator whose action on a smooth function / of the (possibly complex) 
coordinates z = {zi, . . . , z^) is given by 

(Kijf) {Zi, ...,Zi,...,Zj,...,ZN)= fizi, ...,Zj,...,Zi,...,ZN)- (3) 

Given a scalar differential-difference operator D linear in the exchange operator K^j, 
we shall denote by D* the differential operator acting on C°° ® S obtained from D 
by the replacement K^j —>■ Sij. One of the key ingredients in our construction is the 
fact that 

£)$ = £)*$ , for all $ G A(C°° ® S) , (4) 

where A denotes the projector on states totally symmetric under simultaneous per- 
mutations of the coordinates and spins. Consider next the second-order differential- 
difference operators T^ given by 

^-^ ^-^ Zi- Zi+i ^ ^-^ [Zi- Zi+iY 

where di = (9^., z^+i = Zi, and 

'&o{x,y) = l, '&i{x,y) = -{x + y), -dii^x.y) = xy . 

Each Hamiltonian H^ is related to a linear combination 

H, = cT, + c.J- + CoJ'^ + Eo (6) 

of its corresponding operator T^ and the first-order differential operators 

J- = J2^^, J° = 5^z,9, (7) 



TABLE I: Parameters, gauge factor and change of variable in Eqs. ([6]) and 
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through the star mapping, a change of variables and a gauge transformation. More 
precisely, 

H, = ^,■H:\^^^^^^^y^^-\ e = 0,l,2, (8) 

where the constants c, c_, cq, -Eq, the gauge factor /i(x), and the change of variables 
C,{x) for each model are listed in Table [B 

From Eqs. (HI) and (IHl) it follows that if $(z) G A(C°° ® S) is a symmetric 
eigenfunction of if^, then 

x[/(x)=/i(x)<|.(z)|,^=^(.^) (9) 

is a (formal) eigenfunction of H^ with the same eigenvalue. In this paper we shall 
construct a flag 7i° C Til C ■■■ of finite-dimensional subspaces of A(C[z] ® S) 
invariant under each H^. We will show that the problem of diagonalizing H^ in 
each subspace 7i" is equivalent to the computation of the polynomial solutions of 
a system of linear differential equations. We shall completely solve this problem, 
thereby obtaining several infinite families of eigenfunctions of H^ for each e. From 
the expressions for the change of variable and the gauge factor in Table [H and the 
fact that the functions $ in Eq. ([9]) are in all cases polynomials, it immediately 
follows that the eigenfunctions thus obtained are in fact normalizable. 

Remark 2. The operators (I5l) can be expressed as quadratic combinations of the 
first-order operators 

Ui = z^di , Qi = (1 - Ai,i+i) H (1 - Aj-i,j) , 

where e = 0, 1, as follows: 



T2 = j:[iDl)' + {DlaQ] + a-l} 



For each nonnegative integer n, the space P" of polynomials in z of total degree 
at most n is invariant under the operators Ql (see Ref. (35|), and hence also under 

6 



both Tf and H^. Note, however, that the operators H^ do not commute with 
the symmetrizer A, and thus the previous observation does not imply that they 
preserve the space A(P" ® E) of symmetric spin functions of polynomial type. 
Consequently, H^^ is not guaranteed a priori to admit finite-dimensional invariant 
subspaces of A(C[z] ® E) . This is in fact the main difference with the usual solvable 



spin CS models [8|, l30|, l35|, l36|, |4J, |45|, |46[, for which the operators analogous to 



H^ preserve Vn and commute with A, and hence automatically leave invariant the 
space A(P"®E). 



3. THE INVARIANT SPACES 

In this section we shall prove that each operator T^ leaves invariant a flag T^ C 
T^ C • • • , where TJ^ is a finite-dimensional subspace of A('P" ® E). This result 
will then be used to construct a corresponding invariant flag Ti^ C Til C ■ ■ ■ for the 
operator H^, where ?^" C TJ^ for all n. 

Let us first introduce the following two sets of elementary symmetric polynomials: 

^k = ^z^i, rk= Y^ Zi,---Zi^; k = l,...,N. 

i ii<---<ik 

It is well known that any symmetric polynomial in z can be expressed as a polynomial 
in either a = (ai, . . . , a^) or r = (ri, . . . , tj^). 

We shall denote by 2aX^ the terms of T^ linear in derivatives, that is 

^e = y] {ztdi - ^i+i^i+i) . 

"~ Zi — Zi+i 

In the next lemma we show that each vector field X^ leaves invariant a corresponding 
flag X^ (Z X^ G ■ ■ ■ of finite-dimensional subspaces of the space S = C[(t] = C[t] 
of symmetric polynomials in z. 

Lemma 1. For each n = 0,1, . . . , the operator X^ leaves invariant the linear space 
X^, where 

X^ = C[a^, a2, a,] n P^ A'f = C^, a^, r^] n P", X^ = C[a^, t^_^, t^] n P^ 

Proof. If / is a function of the symmetric variables cxi, o"2, a^, t/v_i, t^, we shall use 
from now on the convenient notation 

r 1 - "K^ ' k = 1,2,3, 

' " " k = N-l,N. 

Let us first consider the vector field Xq. Since 




(k—1 
"^-^ Zi- Zi+i "^ 



Zi—\ Z-i 




6o"i , k 



if / G Xq we have 

Xo/ = 2(iV/2 + 3^1/3) G^o"- (10a) 

The proof for the remaining two cases follows from the analogous formulas 

XJ = Nf,+4aj2, feXr, (10b) 

X2/ = 2ai/i + N{TN-ifN~i + tnIn) , feX^. (10c) 

D 

Remark 3. It should be noted that these flags cannot be trivially enlarged, since, 
e.g., 

- XoO"4 = 2(T2 + 2_^ ZiZiJ^i , 

i 

- Xi(J3 = 2(72 + ^ ZiZi+i , XiTn-1 = Tn "^{ZiZi+iy^ , 

i i 

- X2a2 = 2o-2 + ^ ZiZi+i , X2TN-2 = Ntn-2 - Tn ^{ziZi+i)~^ 

i i 

are not symmetric polynomials. 

We note that the restriction of T^ to X^ C S obviously satisfies 

T,\x^^ = J2z'id^ + 2aX,. (11) 

i 

The second-order terms of the operator (llip . however, do not preserve the corre- 
sponding space X^, unless one imposes the additional restrictions specified in the 
following proposition: 

Proposition 1. For each n = 0, 1, ... , the operator T^ leaves invariant the linear 
space 5", where 

'5o" = {/eA'o"|/33 = o}, 

S^ = {fexnf22 = fNN = 0}, 

<S2 = {f ^ "^2 I /ll = fN~l,N-l = 0} . 

Proof. Let us begin with the operator Tq. If / G Xq, an elementary computation 
shows that 

dj = f\ + 2zj2 + ^z}h (12) 

and therefore 

Y, dlf = N{fn + 2/2) + 2(2/12 + 3/3)^1 

i 

+ 2(3/i3 + 2/22)^2 + 12/23^3 + 9/33(14 . (13) 

From the previous formula and Eq. (llOal) it follows that Tq/ G 5o whenever f & Sq. 
Similarly, if / G A"" we have 

dJ = h + 2zj2 + z-W^fN, (14) 



so that 

J2 ^id-f = (/ll + 2/2)ai + 4/12(12 + 4/22^3 

i 

+ 2NfiNTN + 4/2ivairiv + /nnTn-iTn , (15) 

which together with Eq. ( JlObI) iinphes that Ti/ G iS" for all f & S^. Finally, if 
/ G A'a" then 

dif = fl+ {z:f'^TN-l - zr^TN)fN~l + zI'^TnIn (16) 

and hence 

Y. ^'^'f = /ii^2 + 2/i,^_i(air^_i - iVr^v) 

+ 2/iArairAr + fN-l,N-l [{N - l)r^_i - 2r7V-2^Ar] 

+ 2(iV - l)fN^i,NrN-irN + iV/jv7v4 • (17) 
The statement follows again from the previous equation and Eq. fllOcp . D 

The last proposition implies that each operator T^ preserves "trivial" symmetric 
spaces 5" ® A(E) spanned by factorized states. The main theorem of this section 
shows that in fact the latter operator leaves invariant a flag of nontrivial finite- 
dimensional subspaces of A{V^ (S> S). Before stating this theorem we need to make 
a few preliminary definitions. Given a spin state |s) G S, we set 

l^.) = ^E^I^)' l^^^ = M E ^i^±S^,)\s), (18) 

7r(l)=i 7r(l)=i,7r(2)=i 

where S^ is the symmetric group on A^ elements. Here and throughout the paper 
we identify an abstract permutation vr with its realization as a permutation of the 
particles' spins. From Eq. (ITSl) we have 

A{f{z^)\s))=Y,fi^^)\s^), A{g^{z,,Z2)\s))=J29^i^^^^J)\4)^ (19) 

i i<j 

where the last identity holds ii g'^{z2, zi) = ig'^^zi, 2:2). We also define the subspace 

E' = {|s)GS| EI4+i)eA(s)}cs. (20) 

Theorem 1. Let 

V = (/(^i> 0^2, cr3)A|s), 5f((7i, 0-2, a3)A(2;i|s)), /i((7i, (y2)A{z\\s)), 

h{ax,(y'2)k{zxZ2\s')),w{(yi,(y'i)k{zxZi{z\ - Z2)\s)) \ f^s =5-33 = O) , 
'^1" = {ficri,cr2,rN)A\s),g{ai,TN)A{zi\s)) \ /22 = /nn = gwN = O) , 
'^2" = (/(^i7 ^7V-i, i"7v)A|s), 5f(rAr„i, rN)A{zi\s)), r^qiai, Tn)A{z^'^\s)) 

I /ll = fN-l,N-l = gN-l,N-l = Q'll = 0) , 

where \s) G S, |s') G S', deg/ < n, degg < n ~ 1, degh < n — 2, degh < n — 2, 
degw < n — 3, degg < n — N + 1, and deg is the total degree in z. Then TJ^ is 
invariant under T^ for all n = 0,1, ... . 

9 



Proof. By Proposition [T|, it suffices to show that T^ maps T^ / {S"^ ® ^(^)) iiito T^. 
We shall ffist deal with the operator Tq. Consider the states of the form gK{zi\s)), 
with g G Sq~^ . Since 

{di- di+i)zi = {I- Ki^i+i)zi, Wi,l, 

Zl — Zij^i 

we have 

Toigzi) = {Tog)zi + 2dig. 

Calling 

<t>^''^ = A{z'l\s)), keZ, (21) 

from Eqs. ( IT2ll and (IT9l) we obtain 

ro(^$«) = J2Toi9^^)h) = (To^)$« + 2X^A;^fc<|.('=-i) G V'^ (22) 

i k=l 

Similarly, if h{ai,a2) G Sq^'^, the identity 

{di - di+i)zf = (1 - K^+i)zf + {zi - zi+i){5ii + 5i^i-i) , V z, / 

Zl — Zi+i 

implies that 

To{hzl) = {Toh)z^ + AzAh + 2(2a + l)h , 

and therefore 

To{h<S>('^)=J2Uhz^)\s.) 

i 

= {Toh + 8/i2)$^^^ + 4/ii$(i) + 2(2a + l)h¥°^ (23) 

belongs to T^~'^ on account of Eqs. (llOal) and (TT3l) . On the other hand, from the 
equality 

{di - di+i)ziZj = (1 - Ki^i+i)ziZj - {zi - zi+i)5j^i+i5i^i, Mi <j,\/l 

Zl — Zi^i 

it follows that 

To{hziZj) = (TQh)ziZj + 2{zidjh + Zjdih) — 2ah6j^i+i . 

Setting 

¥''> = A{z,Z2\s)) (24) 

and using again Eqs. (TT^ and (I19p we then have 



To{W^)=Y,Uhz.z,)\s±) 

i<j 

= {Toh + Sh)^^^ + 2hiA[{zi + Z2)\s)] ~2ahJ2\4,i+i)- (25) 



10 



Since A[{zi + Z2)\s)] = A[zi{l + ^12) |s)], the RHS of Eq. (^ belongs to Tq^-^ if 
and only if |s) G S'. The last type of states generating the module 7^" are of the 
form w{ai,a2)^^^\ where 

$^^^ = A{zMzi - Z2)\s)) . (26) 

From the equality 

{di - di+i) [ziZj{zi - Zj)] = -^ (1 - i^/,/+i) [ziZ.i{zi - z, 



V-i -1-1-1 y L~i~jv-"i ~jyj / N2 ^^ --1,1 + iy L~«^JV^« ^J^ 

Zl — Zi^i {Zi — Zi^i) 

+ ((5/,i_i + 5ii)zj - {5i,j-i + 5ij)zi, Wi< j,yi 
it follows that 

To[wZiZj{zi - Zj)] = {Tow)ziZj{zi - Zj) + 2zj{2zi - Zj)diW 

- 2zi{2zj - Zi)djW - 2{2a + l){zi - Zj)w . (27) 

Using again Eqs. ( IT2l) and ( TT9l) we obtain 
To(w;8(3)) = Y,To{wz,z,{z, - zj))\sr^) = {T,w + nw^)^''^ 

i<j 

+ 2wiA[izl - zj)\s)] -2{2a+l)wA[{zi-Z2)\s)]. (28) 

Since A[{z^ - 2:3 )k)] = A[zf(l - 5*12)15)], the RHS of the latter equation clearly 
belongs to V'^. This shows that To{%];) C V'^ C Tq". 

Consider next the action of the operator Ti on a state of the form g{ai, tj^)^^^\ 
with g e Si~^. From the identity 

(zidi - zi+idi+i)zi = — (1 - Ki^i+i)zi + -{zi- zi+i){5i^i + (5i,i_i) , Mi, I 

^ Zi — Zi^i z 

we easily obtain 

Tiigzi) = {Tig)zi + 2zidig + 2a5f , 

and therefore, by Eqs. ffTObl) . flTD and (|T5l), 



Ti((7<l>(^)) = Y.T,{gz,)\s,) = {T,g + 2g^)<^^''> + 2{ag + T^gN)<^^'^ G T^' • (29) 

i 

Thus Ti(ri") C r^"-^ C Ti", as claimed. 

Consider, finally, the operator T2. If (^(riv-i, tat) G iSg^^, the identity 

(2;f9i - zf^^di+i)zi = — (1 - Ki^i+i)zi + Zi{zi - zi^x){h,i + (5/,i-i), Vi, / 

zi — Zi+i 

yields 

T2{gzi) = {T2g)zi + 2z,^9i5f + 4azi5f , 

and hence, by Eq. (TTUj) . 

T2((7$«) = 5^T2(<7^,)|.,) 

i 

= [T25 + 2(r^_i(7,v-i + TNgN + 2ag)\ ^^^^ - 2rjv(7jv-i$^°^ (30) 
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clearly belongs to 7^" on account of Eqs. (11 Pel) and (fTTll . The last type of spin states 
we need to study are of the form q^^~^\ where q = Tj\fq{ai, Tj\f) with gn = 0. Since 



{Zi Ul - Zi^^Ul + ijZi = - 

we obtain 



Zi- Zi+i 



) ") 



T2{qzi') = {T2q)zi^ - 2d,q + 2qz-' , 
and thus, by Eqs. flT^ and flTI?|) . 

T2(g<l>(-i)) = ^T2(g^ri)|s,) = (Tag - 2rjvgjv + 2q)<t>^-'^ - 2gi<l'(°) . (31) 

i 

From Eqs. (II Pel) and (TT7|) . it follows that the RHS of the previous equation belongs 
to T^. Hence T2(T^) C 7^", which concludes the proof. D 

Remark 4- We have chosen to allow a certain overlap between the different types 
of states spanning the spaces TJ^. For instance, if \s) is symmetric the state 
g{o'i,a2,(J3)A{zi\s)) G 7^" is also of the form /(o"i, (T2, cr3)A|s). Less trivially, if 
\s) involves only two distinct spin components and is antisymmetric under S'12, then 
we have 

$(3) = |_(^^$(2)_^^$(l))^ 

where ^^^' and ^^^' are respectively defined in Eqs. (I2T1) and (!26l) . Hence, for spin 
1/2 the states of the form w{ai, o"2)$'-^-' in the space 7^" can be expressed in terms 
of the other generators of this space. 

The main result of this section follows easily from the previous theorem: 

Corollary 1. For each e = 0,1,2, the gauge Hamiltonian H^ leaves invariant the 
space 7i" defined by 

Proof. We shall begin by showing that each space T^ is invariant under the operator 
J°. Note first that 

J0$(^)=j$(^), J0$(2) = 2$(2) , jo8(3)^35(3). ^^^^ (33) 

where the states ^^-'^ <l>(^) and <l>(^) are defined in Eqs. (EI]), ((211) and ([26]), respec- 
tively. Using Eqs. (IT2i) . (IT4l) and (IT6l) one can immediately establish the identities 



^V = ^1/1 + 2(T2/2 + 3(T3/3 , V/(ai, a2, aa) , (34a) 

^V = ^1/1 + 2(T2/2 + iVr;v/7V , V/(ai, ^2, r^) , (34b) 

^V = ^1/1 + (iV - l)rN-ifN-i + iVr^/7V , V/((Ti, rjv_i, r^) . (34c) 

From Eqs. (!HH]) - flM]) and the fact that J° is a derivation it follows that J^ leaves 
invariant the spaces T^ for all e = 0,1,2. This implies that H^ preserves T^ for 
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e = 0,2, since the coefficient c_ vanishes in these cases (cf. Table [T]). On the other 
hand, for e = 1 the coefficient c_ is nonzero, and thus we have to consider the action 
of the operator J~ on the space T^. We now have 

j-$(i) ^^-^O-i)^ jGZ, (35) 

and, from Eq. (iT4l) . 

J-f = Nh + 2(Ti/2 + TN-iJN , V/(ai, a2, tn) . (36) 

Hence J~ leaves invariant the subspace 7i" of T^" defined by the restrictions /at = 
gN = 0. From the obvious identity Ti(/<l>(°)) = (Ti/)<I>(°) and Eq. ([29]), together 
with ( JlObI) . ( ITT]) and ( ITSl) . it follows that the operator Ti also preserves 7i". Likewise, 
Eqs. ( I33l) and (134bl) imply that 7i" is invariant under J°, and hence under the gauge 
Hamiltonian Hi. D 

Theorem [T] characterizes the invariant space 7^" in terms of the subspace S' C E 
in Eq. (|20l) that we shall now study in detail. In fact, from the definition of the 
invariant space 7q" it follows that we can consider without loss of generality the 
quotient space S'/~, where \s) ~ \s) ii A{ziZ2\s)) = A{ziZ2\s)). For instance, from 
Eq. ( ITSl) it immediately follows that if |s) G S' and vr G Sn is a permutation such 
that 7r(z) G {1, 2} for i = 1,2, then 7r|s) belongs to S' and is equivalent to \s). 

In the rest of this section, we shall denote |s^) simply as \sij) for the sake of 
conciseness. From Eq. flTSl) it easily follows that any symmetric state belongs to S', 
since 

El^v.+i) = ]^l^)' forall|3)GA(S). (37) 

i 

On the other hand, if \s) G A(S) the corresponding state h{ai,a2)A{ziZ2\s)) is a 
trivial (factorized) state. We shall next show that the reciprocal of this statement 
is also true, up to equivalence. 

Lemma 2. For every \s) G S, A{ziZ2\s)) is a factorized state if and only if \s) ~ 

Ms)- 



iZj 



Proof. Suppose that 

A{ziZ2\s)) = \s)^CijZii 
i<j 

is a factorized state. Since the LHS of the previous formula is symmetric, Cy = c 
for all i,j and |,s) G A(S). By absorbing the constant c into \s) we can take c = 1 
without loss of generality, and therefore 

A{ziZ2\s)) ='^ZiZj\Sij) =T2\s) =^ \Sij) = \s), i,j = 1,... ,N. 

i<j 

From Eq. ( TT9i) with f{zi, Z2) = 1 it then follows that 

A\s) = Y,\s^,) = \N{N-l)\s). 
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Setting I So) = \s) — A|s) and using the previous identity we obtain 

A{ziZ2\so)) = A{ziZ2\s)) - ^.^^^. A|s) = \s)t2 - ^.^^^. A|s) = . 

Hence |s) ~ A|s), as claimed. D 

By the previous observations, it suffices to characterize the nonsymmetric states 
in S'. To this end, let us introduce the linear operator A : S ^ S by 

^|s) = ^|sm+i). (38) 

i 

Given an element |s) = |si . . . s^) of the canonical basis of S, we shall also denote 
by {s^, . . . , s"} the set of distinct components of s = (si, . . . , sat), and by z/j the 
number of times that s* appears among the components of s. For instance, if |s) = 
|— 2,0, 1, —2, 1), then we can take s^ = —2, s^ = 0, s^ = 1, so that z/i = z/3 = 2, 
z/2 = 1- Consider the spin states |xj(s)) = |xj), i = 1, . . . ,n, given by 

\X^) = '^M-^)W's\..) - Y^ Uj{uk-6jk)\s^s''...), z/i>l, (39a) 

l<j,k<Ti 

\X^) = Yl ""3 d^'^' • • • ) + 1^'^' •••))' '^^ = 1 • (39b) 

Here we have adopted the following convention: an ellipsis inside a ket stands for 
an arbitrary ordering of the components in s not indicated explicitly. Note that the 
states ( l39l) are defined only up to equivalence, and that |xj(s)) = Ixj(^s)) for any 
permutation vr G Sat. 

Proposition 2. Given a basic spin state \s) , the associated spin states \Xi{^)) are 
all in S'/~. 

Proof. Consider ffist a state \xi) of the type f l39ap . Using the definition of the 
operator A in Eq. ( l38l) we obtain 

I 'K^Sn — 2 I 

-2E E E ^A^k-5,,)7,\...s^s\..), (40) 

I l<j,fc<?i 7rSS'jv-2 ' 

where the permutations vr act only on the N — 2 spin components specified by the 
ellipses. On the other hand, we have 

+2EE E -A---f---) 



E E E ^.(^.-M^i---f^'---)- (41) 



I l<j,fe<7i 7rGS'jv~2 



14 



Comparing Eqs. (HOl) and (HTl) we obtain 



A\xd = 2(NA\s)-^^J2u,J2 E -l---f---) 






2(a^-2 J] z/j- JA|s) = 2(2z/i-A^)A|s) 

^ l<7<n ^ 



(42) 






This shows that any state of the form fl39al) belongs to S'/~. Suppose next that 
I'i = 1, so that \xi) is given by Eq. (139b|) . Since 



^l^^) = |!EE E ^.^(l•••|V■.••) + l•••f^^••))=4A|s), (43) 



i l<j<n 7reSjv_2 



it follows that in this case \xi) is also in S'/^ 



D 



Remark 5. Just as symmetric spin states, cf. Eq. (|37|) . the states |xj) satisfy the 
relation 

Indeed, if z/j > 1, from Eqs. (I39ap and (l42l) we have 



A|x*) 



z/i(z/i - 1) + J^ 



'^j 






Ms) 



[z/i(z/i - 1) + AT - z/i - (AT - z/i)2] A|s) 
{N-l){2v,-N)^\s) = ^^^A\x^). 



On the other hand, if z/j = 1 Eqs. ( ]39bl) and ( H3l) imply that 

A|x.)=2( Y. ^,)A|s) = 2(Ar-l)A|s) = ^ 



^IXi). 






Example 1. We shall now present all the states of the form ( 139|) for spin 1/2. In this 
case, up to a permutation the basic state |s) is given by 



V N-u 



|s) = ! + ••• + 



')• 



(45) 



If v is either or A^, then n = 1 and thus |xi) is of the type (I39ap and proportional 
to |s). If z/ = 1, then n = 2 and we can take (dropping inessential factors) 



IXi) = 2 (I+- • ■ ■ ) + I-+ ■ • • )) ~ I+- ■ • ■ ) , 1X2) 



)• 



Although the states |xi) and \x2) are linearly independent, the combination 2|xi) + 
(A^ — 2)1x2) is equivalent to a symmetric state. In the case u = N — 1 the states \xi) 



15 



are obtained from the previous ones by flipping the spins. Finally, if2<z/<A^ — 2 
then n = 2 and the states \xi) are now given by 

\Xi) = -\X2) = i^iu-l)\++---)-{N-u){N-u-l)\ ). (46) 



According to the previous example, for spin 1/2 there are exactly n — 1 independent 
states of the form (1391) associated to each basic state |s), up to symmetric states. 
We shall see next that this fact actually holds for arbitrary spin: 

Proposition 3. Let |s) be a basic spin state. If n is the number of distinct com- 
ponents of s, there are exactly n — 1 independent states of the form ( l39l) modulo 
symmetric states. 

Proof. Let p be the number of distinct components s* of s such that z/j > 1. A 
straightforward computation shows that the combination 

n n 

J2\x^)-i2-p)J2^^i^J-^^J)\^'s'■■■)-i^-p)N{N-l)A\s) (47) 

is equivalent to a symmetric state. Suppose first that p 7^ 2. It is immediate to 
check that in this case the set {\Xi) I ^ = 1, ■ ■ ■ ,n} is linearly independent. If a 
linear combination J2i=iCi\xi) is equivalent to a symmetric state \s), then \s) must 
be proportional to A|s), so that we can write 



J2C^\X^)-^K2~P)N{N-1)A\S). 



j=l 



Hence J27=ii^i ~ ^)\Xi) ~ 0, and the linear independence of the states \xi) implies 
that Cj = A for all i. On the other hand, if p = 2 the set {\xi) I i = 1, ■ ■ ■ ,n} is 
linearly dependent on account of Eq. (jlT]), but removing one of the two states with 
I'i > 1 clearly yields a linearly independent set. It is also obvious from the coefficients 
of the states |s*s* . . . ) that no linear combination Yl^=i CjIxj) can be equivalent to a 
nonzero symmetric state. D 



The next natural question to be addressed is whether the states of the form ( l39l) 
span the space S'/~ up to symmetric states: 

Proposition 4. The space (S'/A(S))/~ is spanned by states of the form ( |39l) . 

Proof. For conciseness, we present the proof of this result only for the case M = 1/2. 
Let T,^, denote the subspace of E spanned by basic spin states with u "+" spins, and 
set EJ^ = E' n Ej,. Since the operators A and A involved in the definition (l2Up of E' 
clearly preserve Ej,, it suffices to show that the states |xi(s)) with s given by ( l45l) 
span the space E'j,/~ up to symmetric states. Note first that the statement is trivial 
for i/ = 0, 1, N — 1, N, since in this case the states of the form (l39ll obviously generate 
the whole space Ei,/~ . Suppose, therefore, that 2 < u < N — 2, so that 

E,/~ = (|++...), 1+-...), I )). 
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Since the state (H6ll and the symmetric state (up to equivalence) 

u{u - 1)1++ • • ■ ) + 2z/(iV - z/)|+ ) + (A^ - i^){N - z/ - 1)1 ) 

both belong to S^/~, we need only show that (for instance) \-\ • • ) is not in S^/~, 

i.e., that A\-\ ■ ) is not symmetric. But this is certainly the case, since a state of 

the form 

N—v—k v—k 

k = 1,2, . . . , min(z/, N — v) , 
appears in A\^ — ■ ■ ■ ) with coefficient 2k{v — l)!(iV — v — 1)\ depending on k. D 



4. THE ALGEBRAIC EIGENFUNCTIONS 

In the previous section we have provided a detailed description of the spaces 7i" C 
A(C[z]®S) invariant under the corresponding gauge Hamiltonian H^. In this section 
we shall explicitly compute all the eigenfunctions of the restrictions of the operators 
H^ to their invariant spaces Ti^. This yields several infinite^ families of eigenfunctions 
for each of the models (E]), which is the main result of this paper. We shall use the 
term algebraic to refer to these eigenfunctions and their corresponding energies. It 
is important to observe that the eigenfunctions of the gauge Hamiltonian H^ that 
can be constructed in this way are necessarily invariant under the whole symmetric 
group, in spite of the fact that H^ is symmetric only under cyclic permutations. In 
fact, the explicit solutions of all known CS models with near- neighbors interactions 
(both in the scalar and spin cases) can be factorized as the product of a sirnple 



gauge factor analogous to /x times a completely symmetric function [41|, |42|, |43|, |48 
This, however, does not rule out the existence of other eigenfunctions of the gauge 
Hamiltonian H^ invariant only under the subgroup of cyclic permutations, which is 
indeed an interesting open problem. 

Case a 

We shall begin with the model (l2al) . which is probably the most interesting one 
due to the rich structure of its associated invariant flag. In order to find the algebraic 
energies of the model, note first that one can clearly construct a basis Bq of Tig whose 
elements are homogeneous polynomials in z with coefficients in E. If F G Bq has 
degree k, then J^F = kF and TqF has degree at most k — 2. If Bq is ordered 
according to the degree, the operator Hq is represented in this basis by a triangular 
matrix with diagonal elements Eq + kcQ, where k = 0, . . . ,n is the degree. Thus the 
algebraic energies are the numbers 

Ek = Eo + 2ku, A; = 0, 1, ... . 

We shall next show that the algebraic eigenfunctions of if o can be expressed in 
closed form in terms of generalized Laguerre and Jacobi polynomials. The computa- 
tion is basically a two-step procedure. In the first place, one encodes the eigenvalue 



^ For the model Pc)) we shall see below that the number of eigenfunctions with a given total 
momentum is finite. 
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problem in the invariant space Hq as a system of linear partial differential equations. 
The second step then consists in finding the polynomial solutions of this system. 
Regarding the first step, we shall need the following preliminary lemma: 

Lemma 3. The operator Hq preserves the following suhspaces o/T-Cq: 

7?^,|,^ = (/$(°),^$W,/i$(2)), |s)GS, (48) 

^o,ts)=^o,|s) + P^'^), l^)eS', S,,\s) = \s), (49) 

K,|.) = ^o,|s) + (^$^'^) , l^)es, S,2\s) = -\s), (50) 

where f, g, h, h, w are as in the definition ofT^ in Theorem{l\, and ^^'^^ ^^'^\ ^^^^ 
are respectively given by fl?T]) . ( H^ and 



Proof. The identity To(/<l>(°)) = (To/)<l'(°) and Eqs. ([22]), ([23]) and ([33]) clearly 
imply that the subspace '^^qIs) ^^ invariant under Hq. Consider next the action of 

Hq on a function of the form h^^'^\ Since \s) is symmetric under S12, we can replace 

A [(2:1 + 22) I s)] by 2$^^) in Eq. (^31) . Secondly, any state |s) G S' satisfies the identity 

El<m) = ]^A|.). (51) 

i 

Indeed, we already know that the previous identity holds for symmetric states 
(cf. Eq. ([37])) and for states of the form ([39]) (cf. Eqs. ([38]) and ([B])). On the 
other hand, by Proposition Hj every state in S' is a linear combination of a symmet- 
ric state, states of the form flH^ . and a state \s) such that A(^ziZ2\s)) = 0. But for 
the latter "null" state \sij) = for all i < j, and hence A\s) = A\s) = 0. Therefore, 
Eq. fl^ can be written as 

To (/i$(2)) = (^Toh + 8/12)$^') + 4/ii<l>(^) - -^ /i$(°) . (52) 

From the previous equation and Eq. (!H^ it follows that ifo(^'&*^^-') ^ ^0 |s>- Fiiicilly, 
if 5*1215) = —\s), Eq. f[28]) reduces to 

To(w$(3)) = (TqW + 12^2)$^^^ + 4wi$(^) - 4(2a + l)w<l>(^) , (53) 

which, together with Eq. ([33]), implies that Ho{w$^^^) G 7^^,^^. D 

Remark 6. The requirement that \s) be symmetric (respectively antisymmetric) 
under 6*12 in the definition of the space '^^qIs) (respectively Hqi^J is no real 
restriction, since the antisymmetric (respectively symmetric) part of |s) does not 
contribute to the state $^^'' (respectively $'^^)). 

By the previous lemma, we can consider without loss of generality eigenfunctions 
of Hq of the form 

$ = /$(o) + ^$« + /i$(2) + /i$(2) + ^$(3) ^ deg$ = A;, (54) 
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where the spin functions ^^'^^ $*^^)_and $^^'' are all built from the same spin state 
|s). Note that we can assume that hw = 0, and that the spin state \s) is symmetric 
under Su and belongs to S' if /i 7^ 0, whereas it is antisymmetric under 5*12 ii w ^ 0. 
Using Eqs. ([221), (pD, (ESI), (ESI) and ([53]), it is straightforward to show that the 
eigenvalue equation Hq^ = {Eq + 2ku)^ is equivalent to the system 

[-To + 2cj( J° + 3 - k)] w - 12^2 = , (55a) 

[-To + 2cj( J° + 2-k)]h-8h2 = 0, (55b) 

[-To + 2cj( J° + 2-k)]h-8h2 = 6^3 + Awi , (55c) 

[-To + 2cj( J° + 1 - A;)] ^ - %2 = 4(/ii + h)- 4(2a + l)w , (55d) 

[-To + 2lu{J' - A;)]/ = 2(^71 + (2a + l)h - ^^ /i) . (55e) 

Since / and g are linear in a^, we can write 

/ = P + 0-3? , 5- = M + dat; , (56) 

where p, q, u and v are polynomials in o"i and o"2. Taking into account that the 
action of Tq on scalar symmetric functions is given by the RHS of Eq. (TTTj) with 
e = 0, and using Eqs. OlOal) . (IT^ and 034ap . we finally obtain the following linear 
system of PDEs: 

[To -2uj{k- 3)] w - I2w2 = , (57a) 

[Lo - 2uj{k - 2)]h - 8h2 = , (57b) 

[Lq - 2iu{k - 2)]h - 8h2 = 6v + Awi , (57c) 
[Lq - 2uLi{k -l)]u- 4u2 = 4/ii + Ahi + Qa2Vi 

+ 6{2a + l)aiv-A{2a+l)w, (57d) 

[To - 2iu{k -A)]v- 16t;2 = , (57e) 

4a 
(To - 2ujk)p = 2ui + 2(2a + l)h - — h + 6a2qi + 6(2a + l)aiq , (57f) 

[To - 2ujik -3)]q- I2q2 = 2vi , (57g) 

where 

To = -{Ndl^ + Aa^d^.d^, + 4^29^^ + 2(2a + l)Nd,,) 

+ 2uo{aid„,+2a2d„^). (58) 

As a consequence of the general discussion of the previous Section, the latter system 
is guaranteed to possess polynomials solutions. In fact, these polynomial solutions 
can be expressed in closed form in terms of generalized Laguerre polynomials L^ 
and Jacobi polynomials Pp' . 

Theorem 2. Let 

1^ 3 ..... ,./ 1\ 2r2 



a = N{a + -)--, (3 = (3{m) = l-m-N(a + ^, t 
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where x = j^ Yli ^i ^■^ ^he center of mass coordinate. The Hamiltonian Hq possesses 
the following families of spin eigenf unctions with eigenvalue Eim = Eq + 2a; (2/ + m), 
with I > and m as indicated in each case: 



,(2) 



^1^=/^^'"-%'"^ 



2r-l3(,,,^2-s 



Pg+'f(i)('^^'^-25;<l>«) 



2 p(a+2,/3). 






m > 2, 



-,(2) 



"i-2 r-/3/, ,^2\ 



^1::=/^^'""%""^ 



P:+Y\t)($(2)-2a;$«) 



2(a + l] 



+ x" ( p("+2'/3) (t) + ^^"^'^ p("+l,/3) /^^ \ $(0) 

+ a:lPf^i_, W+(2[^] + l)(iV-l) [fl-i ^V 



m> 2, 



^!2 = f^x-'L;^{ur-) 



$!2 = /.a:--^Lr^(..r^) 



i 

'P^^f\t){m-2x^(^^) 



$(°) , m > 3 , 



+ 2x^ (p^Sf\t) + 1^ P^^;^\t) ] $« 



2[f ] - 1 



' 3 [^] 



w + 



(a+2,/3) 



Pm-3 



2[f]-l [ 2 



(t) 



Mtl = /iX-^V^lcr^) 



m(m — 2) ^- ' 



m > 3 , 



W^] + 2 



i 



m > A. 



Here [■] denotes the integer part, e{m) = (l — (— 1)™')/2, and 

$(*=) = A(x^ I s)), $(2) = A(a;ia;2|s)), $'^^^ = A{xiX2ixi - X2)\s)), 

where the spin state \s) is symmetric under S12 and belongs to S' for the eigenfunc- 

~ (2) ^(3) 

tion "^iJ^ , and is antisymmetric under S12 for the eigenfunction ^i^ . The functions 

^m, (pm and Xm are polynomials given explicitly by 



Vr. 



m + 2a + 2 1^+2,(3-2) p(a+3,/3-i) 4a + 7 (a+2,/3-1) , 1 p(q+3,/3) 



1 ^- 

m — 1 2 



pV«-t-o 

-t m -1 



P^ 

m — I 2 ^ 



- P^ 

3 f -2 ' 
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p(a+4,/3-l) _ ^p(a+3,/3-l) _ TU + 2a + 3 (a+2,/3-1) 

1 m -| '^■^III__1 / \ / \ J^ 211. 1 

2 ^ 2 ^ (m — l)(m — 3) " 



2 



_ 1 p(a+4,/3) m + 2a - 1 (a+3,/3) 

3 2^ m — 3 2 ^ 

3m + 2a ^(a+2,/3-1) 2m — 7 (Q+3,/3-1) p(a+4,/3-i) 



(m — lj(m, — 3j 2 ^ m — 3 2 ^ 



-1 



m + 2a + 2 (a+2,/3) _ ^ + 2a (a+s,/?) 1 p(a+4,/3) 



(m — l)(m — 3) ^~ m — 3 ^" 3 

jor even m, and 

,„ _ 9p{"+2,/3-l) _ p{a+3,/3-l) 1 p(a+3,/3) , ^ + 2a + 2 (a+1,/3) 

V-^m — -^-l m-l -1 m-1 T -T „_3 "T , , 1 rn-3 

~^~ ~^~ 3 ^^ m(?7i — 2) ~^~ 

m + 2a + 2 



P 



("+2,/3) 



2- m — 3 5 

, _ p(a+4,/3-l) _ 2m- 5 p(a+3,/3) _ 1 p(a+4,/3) m + 2a - 1 p(Q+3,/3) 

^^ 771 — 2 ^j— 3 ^j— m — 2 ~^~ 

^ 2m - 3 (a+2,/3-1) 2(m- 3) (a+3,/3-1) _ p(a+4,/3-l) 

m[m — 2) ~^r~ m — 2 ~^r~ ~^r~ 

_ m + 2a + 1 (a+2,/3) _ ^ + 2a (a+3,/3) 1 p{a+4,/3) 

m[m — 2) ~^~ m — 2 ~^r~ 3 ~^r~ 

for odd m. For every n = 0, 1, . . . , t/ie above eigenf unctions with 21 + m < n span 
the whole Ho-invariant space ^Hq. 

Proof. Recall, to begin with, that the algebraic eigenfunctions of Hq are of the form 
\E' = /i$, with fi given in Table [land $ an eigenfunction of -f^o of the form fl3^ - fl3U|) . 
In order to determine $, we must find the most general polynomial solution of the 
linear system (1571) . From the structure of this system it follows that there are seven 
types of independent solutions, characterized by the vanishing of certain subsets of 
the unknown functions p, q, u, v, h, h, w. These types are listed in Table [TTl where 
in the last column we have indicated the eigenfunction of Hq obtained from each 
case. We shall present here in detail the solution of the system ( 157|) for the case 
q = v = h = h = w = and m 7^ 0, which yields the eigenfunctions of the form 
\E'U (the procedure for the other cases is essentially the same). In this case the 
system ( 157|) reduces to 

[Lo ~2uj{k- 1)] u - 4u2 = , (Lo - 2ujk)p = 2ui . (59) 

Let us begin with the homogeneous equation for u. We shall look for polynomial 
solutions of this equation of the form u = Q(o"i, a2)R{cr2), where Q is a homogeneous 
polynomial of degree m — 1 in z and i? is a polynomial of degree / in (T2, so that 
k = deg $ = 2/ + m by Eq. fl5^ . From Eq. fl58l) and the homogeneity of Q we have 

LoiQR) = {LoQ)R + Q{LoR) - 4(TiQii?2 - 8a2Q2R2 
= {LoQ)R + Q(Lo - 4(m - l)d^,)R . 

Hence the equation for u can be written as 

R{Lq - Ad„,)Q = Q{-Lo + 4md^, + 41uj)R, 
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TABLE II: The seven types of solutions of the system (I57p and their corresponding eigen- 
functions. 



Conditions 








Eigenfunction 


q = u = v = h = h = w = 0, 




V 


/o 


Im 


u = v = h = h = w = 0, q 


+ 







Irn 


q = v = h = h = w = 0, u 


/ 







Im 


q = v = h = w = 0, h ^ 








*S 


q = v = h = w = 0, /i/O 








5!S 


q = V = h = 0, w j^ 








Im 


h = w = 0, v^O 








Im 



where Lq = Lo|<^=o- Since (Lq — 4:8^2) Q is a homogeneous polynomial of degree 
m — 3 in z, both sides of the latter equation must vanish separately. We are thus 
led to the following decoupled equations for Q and R: 

(Lo-49,,)Q = 0, (60) 

( - Lo + 4m(9^2 + Alio) R = 0. (61) 

In terms of the variable p = uj(T2, Eq. (1611) can be written as 

4ujC;'^{R) = , 

where 

Ct = pdl + i\ + l- p)dp + V (62) 

is the generalized Laguerre operator. Hence R is proportional to the generalized 
Laguerre polynomial Li{u!a2)- On the other hand, we can write Q = a^~^P{t) 
where P is a polynomial in the homogeneous variable t = "^^p — 1. With this 

substitution, Eq. (lUUI) becomes 

where the Jacobi operator J7i/ is given by 

jM = (1 _ ^2)^2 ^ ^^ _ ^ _ ^^ ^ ^ ^ 2)t] 9i + z/(z/ + 7 + (5 + 1) . 

Thus P{t) is proportional to the Jacobi polynomial P, °„i ' (t), so that we can take 

u = aT-'L;^iuja2)P^^^f\t) . (63) 

We must next find a particular solution of the inhomogeneous equation for p 
in ( l59l) . since the general solution of the corresponding homogeneous equation yields 
an eigenfunction of the simpler type \E'JJ . Since 



u, = <-'^r (^^2) 



n — \ 
2 



{m - l)Pl^f>it) - 2(t + l)P^Zririt) 
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(where the dot denotes derivative with respect to t), we make the ansatz p = 
Q(ci, 0'2)-R(c2), where Q is a homogeneous polynomial of degree m in z and R 
is a polynomial of degree I in a2. Substituting this ansatz into the second equation 
in (l59l) and proceeding as before we immediately obtain 

R{LoQ) + Q{Lo - Amd^2 - AIuj)R = 2ui . 

If we set R = L^ {uja2) the second term of the LHS vanishes, and cancelling the 
common factor L^ {uja2) we are left with the following equation for Q: 



LoQ = 2a 



m~2 
1 



(m - DP, 



ia+1,/3) 



m —- 1 1 
2 J 



{t)-2{t+l)PZ^:^\t) 



The form of the RHS of this equation suggests the ansatz Q = &'^P{t), with P a 
polynomial in the variable t. The previous equation then yields 



r(",/3)/m 



1 



^[f ] ^^> ~ 2N 



[m 



DP, 



(a+1,/3) 



d{"+1./3). 



777 — 1 1 

2 J 



(t)-2(t+i)Pj^ri7(t) 



From the definition of the Jacobi operator we easily obtain 



(64) 



J^ 



{",/3) 



\J \ m — l 1 



+ (l + t)9,--(m-i: 



which implies that P = P °_i is a particular solution of Eq. 



Hence 



P 



)(a+l,/3)/ 



--<Lr^(c.a2)Pj^ri7^(t) 



(65) 



is a particular solution of the inhomogeneous equation in ( 159|) . We have thus shown 
that $ = p$(°) +« <l>(^\ with M and p respectively given by Eqs. ( l63l) and ( l65l) . is an 
eigenfunction of Hq with eigenvalue ii^o + 20; (2/ + m). Multiplying $ by the gauge 
factor jj, we obtain the eigenfunction ^JJ of Hq in the statement. 



It remains to show that the states \E'U {k 



;(2) 



(3) 



0,...,4), ^'J^ and ^ti with 
21 + m < n generate the spaces fHH|) - flHI|) . Consider first the "monomials" of the 
form /icr™cr2$''°\ which belong to yu'^^o u\ if 2/ + m < ra. We can order such mono- 
mials as follows: we say that /iO"™'cr2 $''°^ -< ;ucr™0"2$''°'' if 21' + m' < 21 + m, or 
21' + m' = 21 + m and m' < m. From the expansion [471 . Eq. 8.962.1] 



^i'-'Ht) = ^ E iy (-^);.(7 + 5 + ^ + 1).(7 + A: + 1).-U1 - t)' 



ul ^^ 2^k 

k=0 



where {x)k is the Pochhammer symbol 

(x)fc =x{x + 1) 



(x + k-D 



it follows that Pu' (0) > provided that 7 + 1 > and 'y + 6 + 2u < 0. In particular, 
PtfiO) > since 



a + 1 = N\ a + 



i>iV-i>0, a + /3 + 2 



m 




m 




= 2 




/2 




_y 



1 1 

m < -- < 0. 

2 - 2 
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Hence we can write 

where q^ 7^ 0, so that 

{^f2 \2l + m<n) = (;u<4$(°) \2l + m<n). 
Likewise, a similar argument shows that for m > 1 

Ux'^-^L~'^{ur^)Pl^f\t) $(1) \2l + m<n) = (/xar^V^^^^) | 2/ + m < n) , 

I- 2 I 

and therefore 

(*S> *Im I 2/ + m < n) = (/i<4$(°), /ia™-V^<l>(^) | 2/ + m < n) . 

Proceeding in the same way with the remaining spin eigenfunctions we can finally 
show that 

(vj/J^M A; = 0, . . . , 4 , 2/ + m < n> = /iK^I,^ , 

and that 

/^^o,|s) + (^S|2/ + m<n>=/i7Y^^|,^, |s) G S', S,^\s) = \s), 

f^K,\s) + m2 \2l + m<n)= ^ni\,^ , ^i2|s) = -|s), 

as claimed. D 

Remark 7. By Remark 1, the coefficients of $'•''-' in the spin eigenfunctions \E';J 
and ^;^ yield the two families of eigenfunctions of the scalar reduction Hq^ of the 



model (l2al) presented without proof in our previous paper p|. Earlier work on the 
scalar model H^ had established the existence of two families of eigenfunctions 
of the form ixL^ {ujr'^)p^{^ , with Pi, a homogeneous polynomial of degree i^ > 3, 
only for z/ < 6 and N > v [42]. More recently, Ezung et al. [48|] have rederived a 
very small subset of these scalar eigenfunctions by mapping H^ to N decoupled 
oscillators. 

Case b 

Since Cq = Aoj in this case, reasoning as before we conclude that the algebraic 
energies are the numbers 

Ek = Eq + Akuj, A; = 0, 1, ... , 

where k is the degree in z of the corresponding eigenfunctions of Hi. We shall 
see below that these eigenfunctions can be written in terms of generalized Laguerre 
polynomials. To this end, we begin by identifying certain subspaces of H" invariant 
under Hi. 

Lemma 4. For any given spin state \s) G S, the operator Hi preserves the subspace 

K|s) = (/(^i, ^2) $(°\ ^(^1) $(') I /22 = 0) C 7^^ (66) 

where f and g are polynomials of degrees at most n and n — 1 in z, respectively, and 
^^^^ is given by (^I\i . 
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Proof. The statement follows from the obvious identity Ti(/$''°'') = (Ti/)$'-°'' and 
Eqs. dlOb]), (HI]), (USD, (I2nD, iSD and dSSD- □ 



By the previous lemma we can assume that the eigenfunctions of Hi in T^i i^x are 
of the form 

$ = /$(o) + ^$(1) ^ deg(l> = k<n. (67) 

From Eqs. (1291) . ( l33l) and ( l35l) it easily follows that the eigenvalue equation /Ji$ = 
{Eq + 4A;a;)$ can be cast into the system 



-Ti + cu( J° + 1 - fc) 
-Ti+tu(J°-A:)- (6 



^+2'^" 



^-2^1 = 



1 



J- 



f 



2a 



1 



(68a) 
(68b) 



Since / is linear in a2 (cf. Eq. 



), we can write 
f =p + cTsg, 



where p and q are polynomials in o"i. Using Eqs. fllObp . 
we can easily rewrite the system flUHl) as follows: 



(69) 

(USD, dUbi) and daSD 



[Li-uj{k-l)]g-2gi = 0, 
[Li - u{k - 2)]q - 4qi = , 

(Li -ujk)p= ha + h + -\g + 2Ua + h + -\aiq , 



where 



Li = -ai9^, + 



o-i 



dJO"! 



2a + 6 + -)Ar 



9„ 



(70a) 
(70b) 

(70c) 
(71) 



The last step is to construct the polynomials solutions of the system ( 17CT1) . which 
can be expressed in terms of generalized Laguerre polynomials, according to the 
following theorem. 

Theorem 3. The Hamiltonian Hi possesses the following families of spin eigen- 
functions with eigenvalue E^ = Eq + Akuj: 



^ 



^ 



m 



(0) 
k 

(1) 
k 

(2) 



jjiL 



a-l/, .^2a^(0) 



[ujr 



a+l, 

a+3| 
fc-2' 



ur 



Lor 



k>0, 






k> 1 



X, 



(3r' 



$(0) 



k> 2. 



where a = N{2a + b + l), /3 = A^(4a + 6 + |), and^^^'> = A(xf |s)), with j = 0,1 and 
\s) e S. For eac/i |s) G S anc? n = 0, 1, ... , the above eigenfunctions with k < n 
span the whole Hi-invariant space ^Ti"^,-.. 

Proof. As in the previous case, the algebraic eigenfunctions of Hi are of the form 
\E' = |U<I>, where fi is given in Table [Hand $ is an eigenfunction of Hi of the form flUTl) - 
flU^ . The functions p, g and (/ are polynomials in o"i determined by the system (1 70 p . 
which in terms of the variable t = uai = ur"^ can be written as 

2/3 



^-k-liJ 



ra+3 
^k~2'i 



0, 









9- 



Nu' 



tq . 



(72) 
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where C^ is the generahzed Laguerre operator (cf. Eq. (|62l) ). The general polynomial 
solutions of the first two equations in ( 1721) are respectively given by 

g = c,Llt\{t), q = c,Lltl{t). (73) 

On the other hand, from the elementary identity 

Ct{t'Lltf{t))=l{l + \)t^-'Lltf{t), 

it follows that the general polynomial solution of the third equation in fl72|) is given 
by 

Equations (ff3l) and ( |7^ immediately yield the formulas of the eigenfunctions in 
this case. The last assertion in the statement of the theorem follows from the fact 
that the functions p, q and g in Eqs. (1731) and (1741) are the most general polynomial 
solution of the system (1721) . D 

Remark 8. The spin eigenfunctions ^^^ , j = 0, 1,2, listed in the previous theorem 
are essentially the same as those presented in Ref. [l[ (note that in the latter 
reference there is a typo in the formula for the scalar eigenfunction ipn , namely 
the coefficient a multiplying r^ should be replaced by the parameter (3 defined in 
Theorem E]) . 

Remark 9. It should be noted that for u; = the potentials (j2aj) and (l2b1) scale 
as r~^ under dilations of the coordinates (as is the case for the original Calogero 
model). The standard argument used in the solution of the Calogero model shows 
that there is a basis of eigenfunctions of these models of the form /i(x)L^(a;r^)F(x), 
where F is a homogeneous spin-valued function. The eigenfunctions presented in 
Theorems [2] and are indeed of this form. 

Case c 

The model (Pc|) is of less interest than the previous ones, since we shall see that 
in this case the number of independent algebraic eigenfunctions is essentially finite. 
We shall take, for definiteness, the plus sign in the change of variable listed in Table [T] 
(it will be apparent from the discussion that follows that the minus sign does not 
yield additional solutions). 

Let us first note that the potential for this model is translationally invariant, so 
that the total momentum operator P = —i ^^ d^^. commutes with the Hamiltonian 
if 2- Hence the eigenfunctions of if 2 can be assumed to have well-defined total 
momentum. Equivalently, since 

l^'^-P-H =2J°, (75) 

the eigenfunctions of if 2 can be assumed to be homogeneous in z. Let $ be a 
homogeneous eigenfunction of ii2 of degree k and eigenvalue E, so that \E' = //$ is 
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an eigenfunction of H2 with total momentum 2k (cf. Eq. ( ffSl) ) and energy E. By 
Eq. ( 1751) . the function t^'^ clearly has total momentum 2(/c + NX). In fact, the 
following lemma implies that t^'^ is also an eigenfunction of H2 with a suitably 
boosted energy: 

Lemma 5. Let ^ be a homogeneous eigenfunction of H2 of degree k and eigen- 
value E. Then rj^$ is an eigenfunction of H2 with eigenvalue E + 8k\ + AN\^ . 

Proof. From the identity 

^"^ Zi — Zi+i Z '—^ Zi — Zi+i 

where Di = Zidi, we immediately obtain the following expression for the gauge 
Hamiltonian H2'. 



Zi ~ ^j+1 

ZiZiJ^l 



J {H2 - E,) =Y,Dl + a J] ^i±^ (A - A+i) 



-2-E77^^Tl^(l-^M.i)- (76) 



-An^A 



Since r^ DiT^ = Di + X for any real A, it follows that 

t^^H2Tn = H2 + 8A j° + mx^ . 

Taking into account that J°$ = A;$, we conclude that 

H2 (ri^$) = {E + 8kX + ANX^) (r^^) , 

as claimed. D 

By the previous discussion, in what follows any two eigenfunctions of H2 that dif- 
fer by a power of tn shall be considered equivalent. From Theorem [Hand Corollary [1] 
it easily follows that in this case the number of independent algebraic eigenfunctions 
is finite, up to equivalence. More precisely: 

Lemma 6. Up to equivalence, the algebraic eigenfunctions of H2 can be assumed to 
belong to a space of the form 

n2,is) = {au r^-i, airjv-i, r^) ^^^^ + (1, r^_i) ^^'^ + (1, a,) r^^^-'^ (77) 

for some spin state \s) , where ^'•'^^ is given by ( 1211) . 

Proof. Given a spin state |s), the obvious identity 

T2(/$w) = (r2/)$w (78) 

and Eqs. (E]), ( 130|) . ( 13T1) . and ( l33l) imply that the gauge Hamiltonian H2 preserves 
the space 

?i^,|,^ = (/$(°),^$«,r;vg$(-i^>, (79) 

where /, g and q are as in the definition of 7^" in Theorem [H Let $ G 'H^ig) be 

an eigenfunction of H2, which as explained above can be taken as a homogeneous 
function of z. From the conditions satisfied by the functions /, g and q in (!7U|) and 
the homogeneity of $, it readily follows that $ G t]^?^2,|s) for some /, as claimed. D 
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Theorem 4. The Hamiltonian H2 possesses the following spin eigenf unctions with 
zero momentum 



^3 = Ai 



2a 



2a + - 



_<|,(o) + ^cos(2(xi-Xj))|sj) , ^4 = Ai5Z^^'^(2(^^^~^i))l^i)' 



«7^i 



where |sj) zs defined in (1181) anc? x zs t/ie center of mass coordinate. Their energies 
are respectively given by 

Eo, Ei,2 = ^0 + 4(20 + 1--), ^3,4 = Eo + 8(2a + 1) . 

Any algebraic eigenfunction with well-defined total momentum is equivalent to a 
linear combination of the above eigenf unctions. 

Proof. By Lemma [HI in order to compute the algebraic eigenfunctions of H2 it 
suffices to diagonahze H2 in the spaces 7i2,|s> given by fl77|) . Equations flHUj) . flHTj) . 
(l33l) and ( ITHl) . and the fact that if 2 preserves the degree of homogeneity, imply that 
the following subspaces of 'H2,\s) are invariant under H2: 

(^lr^._l,r;v)$^°^ (80a) 

(ai<|.W), (ai$(°),$«), (airA._i<l>(°),r^$(°),rA._i$«) , (80b) 

(r^_i$(°)), (r^._i<|.W,r;v$(-^)), (air;v-i$(°\ r^v^^^ViT^v^^-'^) • (80c) 

From Eqs. (IHUj) it follows that the alternative change of variables 2;^ = e"^'^* does 
not yield additional eigenfunctions of H2. Indeed, the latter change corresponds to 
the mapping z^ i— > l/z^, which up to equivalence leaves the subspace fl80aj) invariant 
and exchanges each subspace in ( I80bl) with the corresponding one in fl80cp . For this 
reason, we can safely ignore the subspaces fl80cl) in the computation that follows, 
provided that we add to the eigenfunctions of if 2 obtained from the subspaces (180b[) 
their images under the mapping Zk ^-* 1/ z^. 

For the subspaces flHOiD-flSObl). using Eqs. ©, (HI]), (HZD, dSOD, ([33D, and ([McD we 
easily obtain the following eigenfunctions of ii2: 

r^$(°\ E = EQ + m, (81a) 

$('\ E = Eo + 4(2a + l), (81b) 

r^-i$(')-^^^$(°\ E = EQ + A{N + Aa + 2). (81c) 

We have omitted the two additional eigenfunctions 

\ 2a + 1/ 

from the above list, since they are respectively obtained from fl81bp and (I81cp when 
the spin state \s) is symmetric. The eigenfunctions (IHTl) are equivalent to the fol- 
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lowing "zero momentuni" eigenf unctions: 

$^°\ E = Eo, (82a) 

r^^/^$W, E = E, + A{2a + l-^), (82b) 

^^ $(^^ - T^ $(°\ E = i?o + 8(2a+l), (82c) 

where the energies have been computed from those in Eqs. ( IHTl) using Lemma O 
The eigenfunctions of H2 hsted in the statement are readily obtained from these 
spin functions together with the transforms of (182bp and 082cp under the mapping 



Zk ^ l/zk- U 

Remark 10. If the spin state |s) is symmetric, then \si) = ^^'^^ /N for all i, and one 
easily obtains from Theorem|l]the following eigenfunctions of the scalar Hamiltonian 
HI-: 



These formulas agree with those in Refs. [ll] and [48!] (the expression of ips in the 
former reference contains an obvious erratum, while this eigenfunction is missing 
altogether in the latter reference). 



5. SUMMARY AND OUTLOOK 

In this paper we have computed in closed form several infinite families of eigen- 
functions of the spin models with near-neighbors interactions ([2]) introduced in our 
previous paper [l|. Our method is based on the fact that each spin Hamiltonian 
ifg is related to a scalar operator H^ involving difference operators which exchange 
pairs of neighboring particles, cf. Eqs. ([S])-([HD- We have explicitly constructed a flag 
of finite-dimensional polynomial subspaces Ti^ C Ti] C ■ ■ ■ invariant under H^ (see 
Corollary [1]) . For all three models ([2]), we have been able to fully diagonalize the 
gauge Hamiltonian H^ in its invariant spaces 7i". Multiplying each eigenfunction 
of H^ in 7i" by the appropriate gauge factor /x and performing a suitable change 
of variables (cf. Table [T]) we obtain the families of eigenfunctions of H^ mentioned 
above. 

The results obtained in this paper suggest several open problems that we shall now 
briefly discuss. In the first place, it would be natural to study the BCjy counterparts 
of the models ([2]), for which the interaction potential also depends on the sums 
Xi + Xj+i. In fact, in the scalar case this question has already been addressed 



in Ref. 42 . It would also be of interest to construct solvable models with near- 



neighbors interactions of elliptic type, both in the scalar and spin cases; see Refs. [35 
36| for a list of models of this type with long-range interactions. An important 
problem closely related with the subject of this paper is the analysis of the spin 
chains obtained from the models ([2]) by applying the freezing trick. These chains 
are characterized by the fact that the interactions are restricted to nearest neighbors 
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(as in the Heisenberg chain), but their strength depends on the distance between 
the sites (as in chains of Haldane-Shastry type). For this reason, we beheve that 
the study of these new chains could prove of considerable interest. 

Consider, for instance, the chain associated with the model (I2al) . whose Hamil- 
tonian is given by 

i 

where (^i, . . . , ^n) is the unique equilibrium of the scalar potential 
Uo = i:r'^ + y" 7 V7 ^ + T^ 



Z , y^^i •^i—\)\-^i •^%^\) ■ \^i *^2+ly 



2 



in the domain xi < ■ ■ ■ < x^. It can be shown that the chain sites ^j are 
symmetrically distributed around the origin; for instance, for A^ = 4 we have 
^4 = -^^ = (v/3 + l)/2, ^3 = -^2 = (v^ - l)/2. Note, in particular, that the 
chain sites are not equally spaced, as is the case in most spin chains of Haldane- 
Shastry type. In principle, it is not possible to apply the method of Refs. |38l. l39l. l40l 



to evaluate the partition function of the chain (1831) in closed form, since the algebraic 
eigenfunctions of the model (l2al) computed in Section H] do not form a complete set. 
On the other hand, the explicit nature of the algebraic eigenfunctions presented in 
Theorem [2] makes it feasible to compute a number of eigenvalues and eigenfunctions 
of the spin chain flH^ by taking the strong coupling limit a —>■ oo. We emphasize 
that the results thus obtained would be valid for an arbitrary number of spins, and 
thus could be helpful in uncovering general properties of the spectrum. By combin- 
ing this approach with numerical computations for fixed values of N, we expect to 
achieve a reasonable understanding of the properties of this novel type of chains. 
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